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Abstract 

SU{2)l triplet scalars appear in models motivated for the left-right symmetry, neutrino masses 
and dark matter (DM), etc.. If the triplets are the main decay products of the DM particle, and 
carry nonzero lepton numbers, they may decay dominantly into lepton pairs, which can naturally 
explain the current experimental results reported by PAMELA and Fermi-LAT or ATIC. We discuss 
this possibility in an extended left-right symmetric model in which the decay of DM particle is 
induced by tiny soft charge-conjugation (C) violating interactions, and calculate the spectra for 
cosmic-ray positrons, neutrinos and gamma-rays. We show that the DM signals in the flux of 
high energy neutrinos can be significantly enhanced, as the triplets couple to both charged leptons 
and neutrinos with the same strength. In this scenario, the predicted neutrino-induced muon flux 
can be several times larger than the case in which DM particle only directly decays into charged 
leptons. In addition, the charged components of the triplet may give an extra contribution to 
the high energy gamma-rays through internal bremsstrahlung process, which depends on the mass 
hierarchy between the DM particle and the triplet scalars. 
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I. INTRODUCTION 



Recently, the PAMELA satellite experiment has observed an excess in the positron frac- 
tion from 10 to 100 GeV pQ , which confirmed the previous hints from HEAT [2] , CAPRICE 
[3] and AMS-01 [1]. The reports from ATIC balloon experiment showed a rapid rise of 
the total electron and positron flux at a range 300-800 GeV [5] with a peak at around 600 
GeV, which agreed with the PPB-BETS results[6]. More recently, the Fermi-LAT [7J and 
the HESS ^ Cherenkov experiments have released their results for the sum of electron and 
positron flux. Although they did not fully confirm the previous results reported by the 
ATIC experiment, both the measurements indicated excesses over the expected background, 
which suggests the presence of extra sources for the spectra. Besides the astrophysical 
explanations for anomalies by some nearly sources like pulsars and supernova remnants, 
the dark matter (DM) annihilation or decay is one of the most interesting explanation from 
particle physics. 

The PAMELA, ATIC, and Fermi-LAT anomalies may be a consequence of DM particle 
annihilating mainly into lepton final states. However, if the dark matter is generated ther- 
mally, the annihilation cross section obtained from the observed relic density is significantly 
lower than that required by the current data. One has to resort to a large boost factor 
(about 100 — 1000) to explain the observed large positron flux. Note that the most probable 
boost factor from the dumpiness of DM structures is estimated to be less than ~ 10 — 20 
[S]. The large boost factor may come from nonperturbative Sommerfeld enhancement (TU] 
or Breit-Wigner enhancement [HHTS]. An other possibility is that the DM particle may 
slightly decay, and dominantly decay into leptons [H]. In this case the bound from the DM 
relic density is relaxed. The current data require that the lifetime of DM particle should be 
of the order O{10'^^s). 

If the mass of DM particle is very heavy, it may first decay into some lighter intermediate 
states rather than directly into standard model (SM) particles. An interesting intermediate 
state might be SU (2) l triplets, as it is known in many models that they may carry nonzero 
lepton numbers, and do not couple to quarks directly. This leptophilic feature can naturally 
account for the PAMELA results on both the excesses of positrons and the absence of large 
anti-proton flux. 

The stability of the DM particle is usually protected by imposing extra discrete sym- 
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metries. In our previous work [15], we have shown that the fundamental symmetries of 
quantum field theory such as parity (P) and charge-conjugation (C) can be used to stabilize 
the DM particle. We have explicitly demonstrated that in a left-right (LR) symmetry model 
with P and CP only broken spontaneously, a gauge singlet scalar with odd CP parity can 
be automatically stable without imposing any extra discrete symmetries. 

Motivated by the recent experimental results, we consider in this work the possibility of 
DM particle decay by adding soft C— violating interactions into our previous model, which 
will result in the decay of DM particle with SU{2)l triplet scalars as intermediate states. 
The decays of the triplets into the SM gauge boson pairs ly^iy^, W^Z^, Z^Z^ as well as 
the SM-like Higgs bosons h^h^ are all suppressed by the smallness of the vacuum expectation 
value (VEV) of the left-handed triplet as required by the tiny neutrino masses. Thus the 
triplets decaying into quarks indirectly through these states are strongly suppressed. 

If the reported positron excess indeed originates from DM particle decay through the in- 
termediate SU (2) L triplets, the possible associating signals such as the high energy neutrinos 
and gamma-rays are expected. This is because the neutral and singly charged components 
of the triplets couple to neutrinos in the same way as the doubly charged components couple 
to charged leptons. The final state leptons may generate high energy gamma-rays through 
inverse Compton scatterings (ICS) and final state radiations (FSR), which is common to 
many DM models. In this model the doubly and singly charged components of the triplet 
can also produce extra gamma-rays through virtual internal bremsstrahlung (VIB), which 
gives additional contributions to the gamma-ray spectrum at very high energy region. 

In this work we first show how the DM particle decay caused by the soft C-violating 
interactions can naturally explain the observed excess of the positron and electron excesses, 
and then focus on the cosmic-ray signals associated with DM particle decay through in- 
termediate SU{2)i triplets in a LR symmetric model with two Higgs bidoublets [151 116] . 
After exploring the typical parameters which can explain the current PAMELA as well as 
Fermi-LAT or ATIC data, we present predictions for the flux of neutrinos and gamma-rays. 
We find that the predicted neutrino-induced muon flux can be significantly larger than the 
case in which DM particle only directly decays into charged leptons. The energy spectrum 
of diffuse gamma-rays can be enhanced by the VIB processes from internal charged triplets 
as well. Although the analysis is done in a particular model, the conclusions are generally 
valid for other models which involve SU{2) l triplets as intermediate states, such as the DM 
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models in connection to the type II seesaw mechanism for neutrino masses (see, e.g. Ref. 

m)- 

This paper is organized as follows: in Sec. [Tl| we discuss the main features of the model 
and focus on the decay of DM particle induced by the soft C-breaking terms. In Sec. Ill, we 
calculate the DM decay signals which includes the positron fraction, total flux of electron 
and positron, neutrino-induced muon flux and high energy gamma-ray flux. The results are 
summarized in Sec. HVl 



II. DECAYING DARK MATTER IN A LR SYMMETRIC MODEL 

We begin with a brief review of the LR model with two Higgs bidoublets described in 
Ref. p3| [16]. The model is a simple extension to the minimal LR model p77], which is based 
on the gauge group SU{2)l ® SU{2)r ® U{1)b-l- The left- and right-handed fermions 
belong to the SU{2)l and SU{2)r doublets respectively. The Higgs sector contains two 
Higgs bidoublets (2,2*,0), x (2,2*,0), a left(right)-handed Higgs triplet Al{r) (3(1),1(3),2), 
and a gauge singlet 5(0,0,0) with the following flavor contents 




X2 X2 



(1) 



The two triplets Al^r are responsible for breaking the left-right symmetry at high energy 
scale and generating small neutrino masses through the seesaw mechanism. The introduction 
of Higgs bidoublets (p and x accounts for the electroweak symmetry breaking and overcome 
the fine-tuning problem in generating the spontaneous CP violation in the minimal LR 
model. Meanwhile it also relaxes the severe low energy phenomenological constraints. The 
gauge singlet 5* is relevant to the DM candidate. 

The kinematic terms for the scalar fields are given by 

C = TT[{D^<py{Dy)] + Tr[(D^Ai)t(Z}'^Ai)] + (0 ^ L ^ i?) , (2) 

where 

T T 

D^^L = O^Al - tgl^WLf, , Al] - ig'B^AL . (3) 
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Due to constraints from low energy phenomenology, the mixing between Wl and Wr is 
rather small. Thus the SM gauge boson is mostly W^. 

In this work, we shall focus on the scalar sector which is relevant to the DM particle. 
Under the P and CP transformations, the scalar fields transform as follows 





P 


CP 






X 




X* 


(4) 






A* 




S 


S 


S* 





The whole scalar potential contains two parts 

V = Vo + Vi. (5) 

It is required that the dominant part Vq is both P- and CP-invariant while a small part Vi 
contains the soft C-violating interactions. The most general form for Vq is given by 

- Vo = -^fiUS + S*) - filSS* - ^filiS + S*f + V2fi,s{S + S*)SS* 

+-i=^3.(5 + s*f + \s{SS*f - ]ks{s + s*Yss* - ^K{S + s*Y 
6v2 4 lb 



' ^ 1 ~ 1~ 

-^~^^AS + S*) + X^,sSS* - - V(5 + S*f 



i=l 



, (6) 



where 



Oi = Tr(AiAz. + AljA^O,O2 = Tr(0V),O3 = Tr(0t<^ + 0V) 

O4 = Tr(xtx), O5 = 1t{x^x + X^X) ■ (7) 

Due to the C and CP symmetries, Vo only involves combinations of {S + S*) and SS* . 
The interactions containing odd powers of {S — S*) are forbidden as they are P-even but 
C-odd. Furthermore, [S — S*) cannot mix with the Higgs multiplets in Oj because the 
five independent gauge-invariant combinations Oi{i = 1, ... ,5) in Eq. ([T]) are both P- and 
C-even. Other possible Higgs multiplet combinations such as Tr (0^0 — 0^0) and Tr(A|^AL — 
A^A/j) are P-odd, thus cannot couple to {S — S*). The terms proportional to even powers 
of {S — S*) can be rewritten in terms of {S + S*y and SS*. The absence of odd power term 
of {S — S*) means that Sd = {S — S*)/{i\/2) is a potential DM candidate. 



The spontaneous symmetry breaking (SSB) scheme is to reahze the breaking pattern 
SU{2)l ® SU{2)r ® U{1)b^l SU{2)l ® U{1)y ^ f/(l)em. After the SSB, the Higgs 
multiplets obtain nonzero VEVs 

K,n) = = ^ and (x?,,) = ||, (8) 

where Ki^2, ^1,2, "^l and vr are in general complex, and k = + 1^2^ + + 1^2^ ~ 

246 GeV represents the electroweak symmetry breaking scale. The value of vr sets the scale 
of LR symmetry breaking which is directly linked to the right-handed gauge boson masses 
Mwji = gvR/\/2. It is required that Sd does not develop a VEV directly from Vo, namely 
CP is not broken by the singlet sector. It follows that after the SSB, although P and CP 
are both broken, there is a residual Z2 type of discrete symmetry on Sn in the gauge singlet 
sector, which is induced from the original CP symmetry. 

Constraints on the parameter space from the DM relic density has been discussed in Ref. 
|15j . The main DM annihilation channels are SdSd W^W^,tt and SM-like Higgs boson 
h^h'^. Making use of the constrained parameter space, we predicted the weakly interacting 
massive particle-nucleus elastic scattering cross section. The results show that the typical 
scattering cross section is about one order of magnitude below the current experimental 
upper bound. But in the enhanced Yukawa coupling case, the resulting spin-independent 
cross section can reach 10~^'^cm^ which is close to the improved bound set by the recent 
CDMS-II experiment [18]. 

The tiny soft C-breaking terms in Vi may lead to the DM particle decay. The most 
general P-even but C-odd terms with dimension less than four are given by 

" 5 

Y,c^o^ + us + s*y + c7{s-s*) 



-V, = fi,{S - S*] 



1=1 



(9) 



In this case, the Sn may decay dominantly into two lighter scalars such as Sl6l. The lifetime 



of Sd is estimated by 1/td ~ A^e y 1 ~~ 4m2/m|,/(167rm£i) for m <^ m/5, where niD and m 
stand for the mass of Sd and final state scalars respectively. In order to have a lifetime 
around C(lO^^s), the required value of fi^ is around 10"^^ GeV. 

In this work we focus on the case in which (^i ^ Q (i = 2, . . . ,7) and the right-handed 
triplet is much heavier than that of Sd, such that Sd will decay dominantly into left- 
handed triplet Al as shown in Fig. [l| 
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FIG. 1: Feynman diagrams for DM particle Sd decaying into charged leptons and neutrinos through 
intermediate SU{2)l triplets S'l^Sj^^ , S^6j^ and 6^6^. 

The triplets with non-zero lepton number do not couple to quarks directly. The Yukawa 
interaction between the scalars and leptons is 

= - hi{hij(j) + hijcj) + QijX + gijX)^Rj 

- {^iir2^dLj + Ii-iT2/^RlRj) + h.C, (10) 

which leads to the following seesaw formula for left-handed neutrino mass matrix 

M, = Ml-Md^MI (11) 

with {ML{R))ij = yijVL{R)/V2 and (M^,)^^ = (hijKi + hijK* + gij^i + gijQ)/y/2. The upper 
bound for vl from neutrino masses is around a few eV, and that from the p-parameter 
in electroweak precision measurements is about 1 GeV. The absence of Wr from collider 
searches sets a lower bound of vr > (9(TeV). Thus the hierarchies in neutrino masses, SM 
gauge boson masses and the right-handed gauge boson masses require that 

VL<^K<^vn. (12) 

An important consequnce of the smallness of vl is that the couplings between left-handed 
triplets and SM gauge bosons and Z° are strongly suppressed, as df^^W^W^ , dfW^Z^ 
and 5\Z^Z^ couplings are all proportional to vl- For instance, the ratio between the decay 
width of (5++ ^ and 5++ ^ W+W^ is given by 

where Yu is the Yukawa coupling between 5^"*" and charged lepton For the typical 
case Vl ~ 0(10-^) GeV and ms^ ~ 0(10^) GeV, 1(5^+ is always dominant over 

r((5j^^ — )■ W^^W^+) as long as the Yukawa coupling is not too small, namely, Y^^ > 10^^°. 



Similar results are found for the W^Z'^, Z^Z'^ and hb^hb^ final states. Thus a small required 
by the tiny neutrino masses naturally suppresses the triplet non-leptonic decays channels, 
such that our model can avoid the constraints from the PAMELA anti-proton data. 
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FIG. 2: Feynman diagrams for gamma-ray emission through internal bremsstrahlung. (Left) final 
state radiation (FSR) from final state charged leptons. (Right) virtual internal bremsstrahlung 
(VIB) from doubly charged triplet scalars. 

The Yukawa interaction term l£iT2^L^L requires that the couplings Sf^'^V^l^ ^ S^uu and 
S^l^T^ are the same within each generation in the flavor basis. If the positron excesses 
reported by the PAMELA and Fermi-LAT are indeed from the DM decay through triplets, 
there must exist associating high energy neutrino flux which can be detected by the future 
experiments. 

The SU{2) triplet contains doubly as well as singly charged scalars. Besides the ordinary 
ICS and FSR caused by the final state charged leptons, the charged triplet scalars will emit 
gamma-rays through VIB processes in the cascade decay So ^l^^^^I as shown in Fig. 
[2| The photon multiplicity for a decay process S — )■ XX is defined as [19] 

dx CTs^xx dx ' 
where x = lE^j ^ and ^J~s = ms is the center-of-mass energy. The radiation is dominated 
by the collinear photon emission case, which can be approximated by 

dNxx aQ^ fs{l-x)\ 

— F(x)log^^^j (15) 

with F{x) = (1 -|- (1 — xY)/x for fermions and (1 — x)/x for bosons. For the decay process 
Sd (^L^^L^; oiie has s = Tnjj and mx = rns^. For the sequential decay 6^^ — )■ s 
and mx are replaced by m^^ and m£ respectively. For doubly charged particles there is a 
factor of four enhancement relative to that of singly charged ones. The final positron and 
gamma-ray energy spectra of cascade decays, such as — )■ 5^^5]^^ — )■ Ai, depend on the 



mass of the intermediate states, which also provide a possibihty to probe the masses of the 
triplets. 

III. ENERGY SPECTRA OF THE COSMIC-RAY PARTICLES 

The dark matter cascade decay processes So — > Sl6l — )■ 4£, 2£2z/ and 4z/£ contribute to 
new sources of primary positrons, neutrinos and gamma-rays in our galaxy. Besides the mass 
and lifetime of Sd, the final energy spectra of the cosmic-ray particles depend on the mass 
ratio between Sd and the triplets. In this section, we consider two typical mass hierarchies: 
the small mass hierarchy case (SH), in which we take the mass ratio rm = 2ms^/mD — 0.8, 
and the large hierarchy case (LH) in which = 0.1. Note that in the two extreme cases 
with Tm =1 and 0, the shape of the final energy spectra should reduce to that of the two- 
body and four-body decays respectively. The relative strengths of the Yukawa couplings 
also play important roles, as they determine the corresponding branching ratios of different 
decay modes. In the basis in which the left-handed weak gauge interaction is diagonal, the 
leptonic Yukawa couplings Yij {i, j — e, /i, r) are related to i/ij by Y = - y - U, where U is 
the mixing matrix for charged lepton. For the Yukawa couplings, we are going to consider 
four representative cases. In each of the first three cases, the Yukawa couplings are assumed 
to be dominant by one of the three generation leptons, while in the last case the couplings 
are assumed to be the same for all the generations, i.e., I) |Yee| ^ I^mmI' l^r|, the triplet 
scalars decay mainly into 4e, 2e2i/e and Aue', II) ^ |^e|, \ Ytt\, the triplets decay mainly 
into 4/i, 2/x2i/^ and 4i/^; III) \Yrr\ ^ |^ee|, \Y^n\: the triplets decay mainly into 4t, 2T2ur 
and Aut-; and IV) [Yeel ~ l^/t^tl ~ I^ttL the triplet scalars decay to all leptons with the same 
branching ratio. In fact, the case IV will give the same results for the fiavor democratic 
case in which all the Yukawa matrix elements are nearly identical. For simplicity we have 
assumed that the off-diagonal Yukawa couplings Yij {i ^ j) are all small such that no lepton 
flavor violating process are considered. For each case, we consider various values of the mass 
and lifetime of Sd with paying attention to the case that can reproduce the current data. 
Note that in this model, the large mass hierarchy 0.1 is unhkely for the cases LH-I 

and LH-Il unless So is much heavier than 10 TeV, since the lower bounds of the triplet 
masses are expected to be around TcV scale. However, we still keep them for the sake of 
completeness as such light intermediate states may be possible in other models. 
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A. Electrons and positrons 



In the Milky Way, the propagation of cosmic-ray (CR) particles can be approximated by 
a two-zone propagation model, which includes two cylinders centered at the galactic center. 
A cylinder of radius R = 20 kpc with height 2L (L = 1 — 20 kpc) in the z direction delimits 
the CR propagation region. A smaller cylinder with the same radius but with thickness 
2hz = 0.2 kpc models the galactic plane. The solar system is located at Vq = 8.5 kpc 
and Zq = 0. In general, the CR propagation equations contain terms for convection and 
reacceleration effects [20]. However, for electrons and positrons with energy E > 10 GeV, 
these effects can safely be neglected [21] • Thus one only needs to consider the diffusion 
effect. In this case, the steady-state diffusion equation for positron is given by 

- K{E) ■ VVe+ - ^{b{E)U) = Q , (16) 

where fe+{E,r, z) is the positron differential number density in cylinder coordinate and 
K[E) = Kq{E / GeVy is the spatial diffusion coefficient. Through fitting the measured ratio 
of boron to carbon (B/C), the propagation parameters 6, Kq and L can be determined 
|22j . In Table. |l| we list three typical combinations of the propagation parameters [211 123]- 
The parameter set MIN and MAX corresponds to the minimal and maximal positron flux 
respectively, while the MED scenario best fits the B/C ratio and produce the moderate flux. 
The energy loss rate b{E) are mainly due to the ICS and synchrotron processes. The two 
processes combined give b{E) = i?^/(GeV ■ te) with te = 10^^ s [21]. For decaying DM 
model, the source term Q is given by 

Q(r,E) = ^j:T/^, (17) 

niD dE 

k 

where Tk is the decay width and k runs over all the channels with positrons in the final 
states. dn^+/dE is the positron energy spectrum per DM decay. The DM halo profile p(r) 
is usually parameterized as a spherically symmetric form 

with local DM density pQ = 0.3GeVcm~^. Possible DM halo profile parameters a,/9,7 and 
Ts are listed in Table. |TTj In this paper, we adopt the NFW profile and MED propagation 
model in numerical calculations. 
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Propagation 5 Kq (kpc^/Myr) L (kpc) 



MED 0.70 0.0112 4 

MAX 0.46 0.0765 15 

MIN 0.55 0.00595 1 



TABLE I: Three typical combinations for the propagation parameters 6, Kq and L from the 
measured B/C ratio [211 [23]. 

Halo Profile a /3 7 (kpc) 
NEW 13 1 20 

Isothermal 2 2 5 
Moore 1.5 3 1.3 30 



TABLE XL The values of parameters a,/3,7 and r<j for the NEW [25], Isothermal [26] and Moore 
[27] DM halo profiles. 



Eq. (16) can be solved analytically from the above considerations. For the decaying dark 
matter, the solution of the positron diffusion can be written as [23] [28] 



1 rmo 

U{E,r,z) = — / dE'Ge+{E,E',r,z)y2T, 
mo Je V 



dE' 



(19) 



The explicit form of the Green's function Ge+{E, E' ,r, z) is given by 



G,+ {E,E',r-z) = J2 G",™^Jo(c 



n,m=l 



xexp 



KqTe 



r \ 


~mTC 








.2L ^ 






4L2 ; 





with 



Gr 



R 



dr'r'Jo Cn 



R 



[z - L) 
5-1 

dz p{r' , z) sm. — —{z' — L) 
L l2L 



(20) 



(21) 



where Jo(i) is the zeroth(first)-order Bessel function, (n are successive zeros of the function 
Jo- At the heliosphere boundary (r = and z = 0), the interstellar positron flux per unit 
energy from the DM decay is 



(22) 
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The interstellar e and background is approximated by 



1 + ll£;o-9 + 3.2^215 
0.7^0-7 



1 + 110^1-5 + 600^2-9 + 580^4-2 
4.5^0-7 



1 + 650^2.3 + 1500^4-2 ' ^^^^ 

where E and $ are in units of GeV and GeV~^cm~2s^^sr^^, respectively. The interstellar 
(IS) positron fraction is given by 

where we have introduced two real parameters /c+ and k_ to normalize the positron and 
electron background, which reflects the background uncertainties [29]. In the numerical 
calculations we take k+{k ) = 0.9(0.7). 

For the solar modulation effects we adopt the Gleeson and Axford analytical force-field 
approximation in which the final flux $(-E©) at the top of atmosphere is given by 

2 

<l>(Ee) = ^HEis) , (25) 
Pis 

where ^{Eis) denotes the primordial interstellar flux. Ejs and E^ have the following relation 

Ejs = E^ + \Ze\(t)F , (26) 

where \Z\ = 1 is the charge number. For the solar modulation parameter (pp, we take 
a typical value of ^i;' = 500 MV. We assume that the solar modulation effect is the same 
for electron and positron, which leads to PF{E^) = PF{Ejs)- Note that even in the case 
of charge-independent modulation, the solar modulation effect leads to a translation for the 
positron factor PF, thus must be considered. 

Making use of the above mentioned formulas we calculate the positron fraction and the 
total electron and positron flux for the eight cases SH(LH)-I~IV. The injection spectra 
dn^/dE is obtained by using the package Pythia 8.1 ^U\. In each case we consider some 
typical values of the parameters m/? and td which can reproduce the experimental data 



well. The favored values for the parameters are listed in Tab. |III| and the corresponding 
numerical results are given in Fig. [3] From the case SH(LII)-I to SII(LII)-III, the required 
mass of Sd increases from 2 to 8 TeV while the lifetime decreases from r/j = 1.5 x 10^® s to 
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case 



ruD (TeV) td (lO^^ s) 2ms^/mD 



SH(LH)-I 2.0 

SH(LH)-II 4.0 

SH(LH)-III 8.0 

SH(LH)-IV 2.5 



1.5 0.8 (0.1) 

0.9 0.8 (0.1) 

0.4 0.8 (0.1) 

1.3 0.8 (0.1) 



TABLE III: Typical values of the mass and lifetime of So as well as favored by the PAMELA 
and Fermi-LAT (or ATIC) data in different cases from SH(LH)-I to SH(LH)-IV. 




FIG. 3: The predicted positron fraction (left) and the total electron and positron flux (right) for 
four representative cases from SH-I to SH-IV. The relevant parameters are given in the figures in 
compact notations. For instance, the notation "SH-II-4.0-0.9" corresponds to the case of SH-II 
with mD=4.0 TeV and td = 0.9 x lO^^s. 

To = 0.4 X 10^^ s. The case SH(LH) -IV is a kind of mixture with equal Yukawa couplings 
for all generation leptons, and the required mass and lifetime are found close to the case 
SH(LH)-I. As indicated in Fig. |3| the SH-I case can explain both the PAMELA and ATIC 
data. In this case, the DM particle decays directly into 4e, which usually leads to a too hard 
spectrum to meet the Fermi-LAT data. The SH-II~IV cases involving fi- and r-lepton as 
final states can reproduce the PAMELA and Fermi-LAT results well with the parameters 



in Tab. |III[ The differences among the four cases are not obvious in the low energy rang 
E < 100 GeV. However, in higher energy range E >! TeV, the heavier final state case leads 
to harder energy spectrum. The SH-III case with 4r final states has the hardest spectrum 
among all the cases in this energy region. 
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FIG. 4: Effects of the DM particle mass and lifetime on the positron fraction (left) and the total 
electron and positron flux (right) for four representative cases from SH-I to SH-IV. For each case 
the mass and life-time of So are varied by a factor of two. 



The mass and hfetime of Sd are two key parameters in determining the positron fraction 
and the total electron and positron flux. In Fig. [4| we vary the DM mass or lifetime by a 



factor of two for each case listed in Tab. Ill The results show that in general the location of 
the peak in the spectrum is controlled by m/^, and the height of the peak is more relevant to 
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FIG. 5: The same as Fig. [4| but for different mass and lifetime of So which can have nearly 
the same results for the positron fraction but significantly differ in the total flux of electron and 
positron. See text for detailed explanation. 



tb. Compared with the PAMELA data, the Fermi-LAT data is more sensitive to mo^ which 
is due to more data points at high energy region and higher statistics. If simultaneously 
adjusting both rriD and r^^ in all the cases from SH-I to SH-IV, one can obtain the nearly 
the same positron fraction [E ^ 100 GeV) to explain the PAMELA data as shown in Fig. 
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FIG. 6: Difference between tlie small and large hierarchy scenarios in the positron fraction (left) 
and the total flux of electron and positron (right) for the cases SH(LH)-I and SH(LH)-II. 

[5j However, the predicted total electron and positron flux are obvious different. 

For a comparison between SH and LH cases, we show in Fig. |6]the results for SH(LH)-I 
and SH(LH)-III. It is clear that the differences are rather mild for the parameters given 
in Tab. |III[ Actually, the SH and LH cases have significant differences in the injection 
energy spectrum dn^+/dE. However, the differences are smeared out due to the energy loss 
processes from the ICS and the synchrotron. 



B. Neutrinos 



The expected muon neutrino flux coming from the galactic center (GC) and arriving at 
the Earth can be estimated by [31] 



Pere 



dn^ 



dE„ 



\a=e,iJb,T k 

where the neutrino oscillation probabilities P^^^u^ are Pu^^Uf, = 0.22, P^^ 



(27) 



0.39 and 



P. 



0.39 [32]. One can have the same results for the muon anti-neutrino flux $f7„. The 



averaged number flux Jah in a cone with half-angle around GC is given by 

1 1 



J. 



An 



dn 



p{l)dl 



(28) 



P&^Q A^ JaO J LOS 

where Af2 = 27r(l — cos0) is the solid angle. The above equation can be written as the 
following form 

1 27r 



J. 



AQ 



PqTq AQ 



d cos d)' 



P 







2/r„ cos0' + /2 dl 



(29) 
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FIG. 7: The predicted neutrino-induced upgoing muon flux as a function of the cone half- angle 
around the galactic center, (left) the direct and secondary neutrino contributions in the SH-II case. 
(Right) the predicted total neutrino-induced upgoing muon flux for the cases SH-I ~ SH-IV. 



MW 



Tq sin^ (f)' + r(T) cos 0' and the DM halo 



where the integration upper limit / 
size tmw ~ 30 kpc. 

The muon neutrinos interact with the earth rock to produce the upgoing muon flux, 
which can be detected by the water Cherenkov detector Super-Kamiokande (SK) |33l 
The neutrino induced muon flux is give by [35] 

ria 2^ 

a=p,n 



(30) 



where L{En) is the range of a muon with energy E^^ until its energy drops below the SK 
threshold Ethr = 1.6 GeV: 



HE,) 



(31) 



where a,. 



p(3^, a, + f3f,Ethr 

2 X lO^^g"^ GeVcm^ and = 4.2 x lO^^g^^cm^. p is the density of the 



material in gem ^ and ^ rap/rrip are the number densities of neutrons and protons with 



5/9 and r„ ^ 4/9. For the detection cross section, we use 
rf<(E.J ^2m,G^ (^a^j^a El 



dE„ 



TT 



'E?. 



(32) 



where Ay = 0.25,0.15, By = 0.06,0.04 and 

As mentioned in Sec II, an important feature of this model is that the Yukawa couplings 
~, (5j£~z/^ and (5°i/^z/£ are equal. In the DM particle decay processes, the final state 



for (5++ 
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neutrinos may come from the direct decay of Sl and from the secondary decay of charged 
leptons. If a large decay rate of 6^^ — t- is needed for explaining the current experimental 
data, the associated processes — )■ i'^i'i and 5° — > uiu^ will give extra contributions to the 
neutrino flux, which can be more important than that neutrinos from the secondary decay 
of charged leptons. In Fig. [7| we compare the contributions to the total neutrino-induced 
muon flux from the direct neutrinos and that from the secondary lepton decays in the SH-II 
case. The difference can be clearly seen in the large angle region (p ^ 30° in which the direct 
neutrino contribution is roughly five times as many as that from secondary charged lepton 
decays. In Fig. [7| we also give the predicted muon flux for the cases SH-I ~ SH-IV. The 
results show that the SH-III case with Sd decaying into 4r, 2T2ur and predicts much 
larger flux than the other cases. The reason is that the direct neutrinos in this case have 
higher energy than that in other cases. This case can be tested in the future experiments 
with improved sensitivity by a factor of a few. Note that the 2r or 4r final state cases are 
now strongly disfavored by the latest Fermi-LAT gamma-ray data [M] even in the DM decay 
scenario. The predictions for the cases LH-I ~ LH-IV are found to be very similar. 



C. Gamma-rays 

The high energy gamma-rays generated by DM particle may come from two type of 
contributions. The first one is that the DM particle directly produces high energy gamma- 
rays through processes VIB and FSR. The second contribution comes from the ICS process. 
In this case, the low energy interstellar photons obtain energy from the high energy positrons 
and electrons. For the first contribution, the expected spectrum of high energy gamma-rays 
is given by 



dE^ 4:71171 D 




^AJVr,-! , (33) 



where Jao is defined in the same way as for neutrino in Eq. (28). For a general observable 
region, it is convenient to choose the galactic coordinate (6, /) where h and / stand for the 
galactic latitude and longitude. Then J^n can be written as the following form 

J An = — — -Jt^ [ dQ [ p ( \/rl. - 2srQ cos 6 cos / + ) ds , (34) 
PeTQAQJ^^ Jo \^ J 

where the integration upper limit Smax = \/"''\iw ~ '^©(l ~ ^^^^ fecos^ /) + cos h cos / and 

dVL = cos bdbdl. 
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The high energy gamma-rays produced by the ICS process e^7 — )■ have the following 
energy spectrum [371 EH] : 

^ il/^L/V/^-'^"^-^' ^'-tt ■ '''' 

Since the photon may be emitted from both electron and positron, an overall factor of 2 has 
been multiplied in the equation. The differential energy density u^{E^,r, z) of interstellar 
radiation field (ISRF) contains three components: the cosmic microwave background (CMB), 
thermal dust radiation and star light. Here we adopt the GALPROP numerical result for 
in Ref. [39]. The parameter fjcs is defined by [37] 

fjcs = 2glng + (1 + 2g)(l - q) + ^^^(1 - q) , (36) 



with 



Ej' ^ ^ Eyml ^^^^ 

Ee AE^Ee{Ee — Ey) 



where < g < 1. In Eq. (35), fe+{Ee,r,z) is the positron differential number density. 



Because of the ICS and synchrotron, high energy electrons and positrons will loss most of 



their energy within a kpc. Therefore one may neglect the diffusion term of Eq. (16) and 
approximately calculate fe+{Ee,r,z) for every point in the CR propagation region. In this 
case, fe+{Ee, r, z) is given by the following formulas [371 SO] 

bics{Ee,r,z) mo ^ Je, dE' 
where the electron energy loss rate bjcs{Ee,r, z) is 

7 

Here we neglect the synchrotron energy loss rate bgyn as bsyn ^ bjcs [3H]- It is worthwhile 
to stress that the above energy loss rate bjcs{Ee, r, z) is position dependent ^ 

In order to compare with the experimental data, one needs to know the diffuse gamma- 
ray background which includes a galactic <|)Gaiactic contribution and an extragalactic (EG) 
contribution. The galactic gamma-ray background <|)Gaiactic j-Qaj^iy comes from pion 



b,cs{Ee, r, z) = j rfE /^^^;^'^^ j dEy{Ey - E,)fjcs . (39) 



^ If a position independent energy loss rate is assumed, one can use Eq. (191 to calculate f^+{Ee,r, z) 
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FIG. 8: The predicted gamma-ray spectra in the small hierarchy (left) and large hierarchy (right) 
scenarios for the observed region 10° < |6| < 20°. 

decay, ICS and bremsstrahlung. In principle, one should run the GALPROP code [20] 
to obtain <|)Gaiactic given diffusion parameters S, Kq and L. Note that <|)Gaiactic jg 
sensitive to these given diffusion parameters as shown in Refs. [H]. For an illustration, we 
use the numerical results of the conventional GALPROP model (44_500180) in Ref. [12] as 
an estimate of our galactic gamma-ray background <|)Gaiactic_ p^^, ^j^g diffuse extragalactic 
gamma-ray background, one can adopt the following parametrization from the Fermi-LAT 
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FIG. 9: The predicted gamma-ray spectra for the observed regions 0.25° < \h\ < 2.75°, 0.25° < 
|/| < 2.75°, 357.25° < |6| < 359.75° (left) and 60° < |6| < 90°, 0° < |/| < 360° (right) in the small 
hierarchy scenario for the cases SH-I ~ SH-IV. 



preliminary results 



GeV J 



-2.45 



(40) 



In Fig. |8| we show the contributions from ICS, FSR and VIB in various cases in the sky 
region 10° < |6| < 20°, 0° < / < 360°. For all the cases in lower energy region E < 100 GeV, 
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the DM contributions are dominated by the ICS process which is typically 1 ~ 2 order of 
magnitude below the current Fermi-LAT data [43j . The ICS process contributes to a broad 
spectrum from 10^^ GeV to a few hundred GeV. At higher energy region ~ 10 — 100 
GeV, the ICS contribution becomes significant as the background drops rapidly. The ICS 
has no significant dependence on the mass hierarchy as it only depends on the spectrum of 
the final state electrons. In general, compared with ICS processes, VIB and FSR contribute 
to gamma-rays at higher energy. For photon energy above 500 GeV, they become dominant 
sources and may lead to an up turn of the photon spectrum. This prediction can be tested 
by the future gamma-ray detectors. The similar conclusions can be obtained for the galactic 
center region (0.25° < |6| < 2.75°, 0.25° < / < 2.75°, 357.25° < / < 359.75°) and the galactic 
pole region (60° < \b\ < 90°, 0° < |/| < 360°). As shown in Figs, gandg the DM particle 
decay can give significant contributions to the high energy gamma-rays for the three typical 
regions. For an illustration purpose, we naively sum up the DM contributions (ICS, FSR 
and VIB) and the diffuse gamma-ray background ($Gaiactic ^^'^). One finds that all the 
cases are compatible with the current Fermi-LAT preliminary results 
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FIG. 10: Comparison of the predicted energy spectra between FSR and VIB process. The curves 
corresponds to Sd = 2.5, 4.0 and 6.0 TeV respectively, with the mass of triplet scalar ms^ fixed at 
1 TeV and td = 0.9 x lO^^s. 



The FSR and VIB processes show stronger dependence on the masses of the triplet scalars. 



and the two contributions are correlated. It follows from Eq. (15) that for fixed masses of 



initial DM particle and final states, the decrease of the triplet mass will enhance VIB while 
shghtly suppress FSR. Thus the VIB process can be important for the case with large mass 
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hierarchy between Sd and Sl- In Fig. 10, we give the FSR and VIB contributions to the 
gamma-ray spectrum for three different DM mass cases (m^i = 2.5,4 and 6 TeV). Here we 
have fixed the triplet mass m^^ = 1 TeV and the DM hfetime td = 0.9 x 10^^ s. We can see 
that with rriD increasing, the VIB contribution becomes more significant relative to FSR. 
For a low mo = 2.5 TeV, the contribution from VIB relative to that of FSR is only ~ 7%. 
But for a larger niD = Q TeV, the contribution can reach ~ 20% which is non-negligible. 



IV. CONCLUSIONS 



We have discussed the DM cascade decay induced by tiny soft C-violating interactions 
in an extension of a left-right symmetric model in which the DM particle is identified as a 
gauge singlet scalar. In this scenario, the DM particle may decay favorably into leptonic 
final states through the intermediate SU {2)l triplets. We have explored the parameter space 
which can well explain the current PAMELA and Fermi-LAT or ATIC experimental data. It 
is shown that this scenario predicts significant associating signals in the flux of high energy 
neutrinos and gamma-rays, which is due to the correlated couplings to neutrinos and the 
extra contributions from the charged components of the triplet. We have found that the 
predicted neutrino-induced muon flux is dominated by the contribution from the neutrinos 
directly from the triplet decay. The gamma-ray radiation from the charged triplet scalars 
through VIB process can be significant in the case that the triplets are much lighter that 
the DM particle. 
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